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Abstract. The usefulness of Riemann P-symbols in deriving identities involv- 
ing the parametrized special function HI is explored. HI is the analytic local 
solution of the Heun equation, the canonical second-order differential equation 
on the Riemann sphere with four regular singular points. The identities dis- 
cussed include ones coming from Mobius automorphisms and F-homotopies, 
and also quadratic and biquadratic transformations. The case when HI is 
identical to a generalized hypergeometric function of 3F2 type is examined, 
and Pfaff and Euler transformations of 3^2(11 , «2, e + l;&i,e;a;) are derived. 
They extend several 3F2 identities of Bailey and Slater. 



1. Introduction 

The Gauss hypergeometric equation (GHE) and the Heun equation (HE) 
are canonical second-order Fuchsian differential equations on the Riemann sphere 
with three and four singular points, respectively. (A differential equation is Fuchs- 
ian if each of its singular points is regular.) By convention, they are written as 
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respectively. Here a,/3;7, S, resp. a(^ 0, 1), q; a, (3; 7, S, e are complex-valued pa- 
rameters; and each equation is invariant under a <-» (3. The singular points are 
x = 0,l,oo, resp. x — 0,1, a, 00, and their characteristic exponents are 0,1 — 7; 
0, 1 — 8; a, (3, resp. 0, 1 — 7; 0, 1 — 5; 0, 1 — e; a, j3. These parameters are constrained 
by Fuchs's relation on exponents. The sum of the exponents of the GHE, resp. 
the HE, must be 1, resp. 2. So 5 — a + [3 — 7 + 1, resp. e — a + j3 — 7 — 5+1. If 
this Fuchsian condition does not hold, the singular point x = 00 will be irregular. 
The parameter q G C of the HE is 'accessory': it does not affect the exponents. 

If 7 is a nonpositive integer, the local (Frobenius) solution of the GHE or HE 
corresponding to the zero exponent at x = will gcncrically be logarithmic; we do 
not consider this case. If 7 is not a nonpositive integer, the solution will be analytic 
at x = 0. When normalized to unity at a: = 0, it is denoted 2F1 (x) := 2-^1(0;, /3; 7; x), 
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resp. Hl(x) := Hl(a,q;a, /?; 7 ,S;x), and is called the Gauss hypergeometric func- 
tion, resp. the local Heun function. (The shorter term 'Heun function' has a more 
specialized meaning |18j .) The redundant parameter S, resp. e, is suppressed. 

The parametrized special functions 2F1, HI have series expansions of the form 
c(n)x n . The coefficient sequences c: N — > C satisfy respective recurrences 

(2a) |(n + j)(n + l)E - (n + a)(n + /?)} c = 0, 

(2b) |(n + 7 + l)(n + 2)aE 2 - [(n + l)(n + 7 + S)a + (n + l)(n + 7 + e) + q] E 

+ (n + a)(n + /3)Xc=0, 

where E denotes the shift operator, i.e., (E k c)(n) := c(n + k). The recurrence (l2"al) . 
resp. (f2b|) . is initialized by c(0) = 1, resp. c(0) = 1, c(— 1) = 0. If c(n) is taken to 
equal zero when n < 0, then (j2aj) . (|2b)) will hold for all n E Z. The power series 
converges on |x| < 1, resp. |x| < min(l, \a\). 

The special functions 2F\,Hl play a central role not only in the theory of 
second-order Fuchsian differential equations, but also in the theory of linear recur- 
rences with quadratic coefficients. From (|2"a)) . (|2b[) , one easily deduces the following. 

Proposition 1.1. 

(1) Suppose c: N — > C (with c(n) := for n < by convention) satisfies a 2- 
term recurrence [P\(n)E + Po(^)] c = 0, where degPi = degPo = 2, and 
one root of Pi is —1 and the other is not an integer greater than — 1 . Then 
the ordinary generating function (o.g.f.) of c equals C 2F\(a, /?; 7; Ax) for 
some A(^= 0), C G C and some choice of GHE parameters a, 0; 7. 

(2) Suppose c: N — » C (iot£/i c(n) := /or n < by convention) satisfies 
a 3-term recurrence [P2(n)E 2 + P\(n)E + Po(n) \ c = 0, where degP2 = 
degPo = 2, degPi ^ 2, and one root 0/P2 is —2 and the other is not an 
integer greater than —2. Suppose, moreover, that the characteristic poly- 
nomial of this recurrence has distinct roots. ( This polynomial is formed 
from the coefficients of n 2 in P 2 , Pi, Po-' it equals p 2 .2n 2 + P\ : 2 n + Ps,2-) 
Then the o.g.f. of c equals C Hl(a, q; a, [3; 7, 5; Ax) for some A(^= 0), C S C 
and some choice of HE parameters a, q; a, [3; 7, S. 

In part (2), the taking of the polynomial roots to be distinct avoids the case when 
the differential equation satisfied by the o.g.f. is confluent, i.e., non- Fuchsian. 

We explain here how the classical formalism of Riemann P-symbols (see §[2]) 
facilitates the derivation of Heun identities, such as alternative expressions for the 
local function HI, written in terms of itself. In §[3J we begin by determining the 
automorphism groups of the GHE and HE. We showed in a previous treatment |14j 
that the two automorphism groups, of orders 24 and 192 respectively, are the Cox- 
eter groups D 3 (= S4) and D 4 . (The set of 24 GHE automorphisms was first 
worked out by Kummer, but its group structure long remained obscure; an isomor- 
phism to the octahedral group S4 was pointed out by Dwork and others [5lll2j.) 
D3 and D4 have subgroups D2 and D3, which are the transformation groups of 2P1 
and HI. The former, of order 4, is generated by classical transformations of Pfaff 
and Euler, but the latter, of order 24, is a novel object. In [14) . we computed the 
24 resulting expressions for Hl(x) with the aid of a computer algebra system. We 
show here how with the aid of P-symbols, they may be computed by hand. 
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In §01 we compute quadratic and quartic transformations of HI, which express 
Hl{x) in terms of Hl(R(x)), where R is rational of degree 2 or 4. It is known that 
the derivation of the quadratic transformations of 2F1 is facilitated by the P-symbol 
formalism 2, §3.9]. Erdelyi [8] §15.3] pointed out that it can also be used to derive 
quadratic and quartic HE transformations, but did not compute explicit formulae 
involving HI. Up to automorphism, we find one quadratic transformation of HI 
and one quartic one that is biquadratic, i.e., is the composition of two quadrat- 
ics. The HI situation apparently contrasts with that of 2-F1: not all the quartic 
transformations of 2-F1, which were worked out by Goursat [9], are biquadratic. 

Using P-symbols, we investigate in §[5] the interesting case when HI coincides 
with a generalized hypergeometric function, in particular with a 3F2 of the form 
3^2(01, a2, e + 1; bi, e; x). If the parameters a, q; a, f3; 7, S are suitably chosen, this 
phenomenon will occur, as noted by Letessier, Valent, and Wimp [llj . Since the 
coefficients of the series expansion of 3F2 satisfy a first-order recurrence, the equal- 
ity is attributable to a factorization of the second-order difference operator in (|2b[) . 
Such factorizations can be computed algorithmically [6], 116] . We give an explicit 
one, and in the same spirit, factor the third-order differential operator in the equa- 
tion satisfied by the 3-F2. The operator products in the two factorizations can be 
viewed as desingularizations in the sense of Abramov, Barkatou, and van Hoeij [1] . 

In §[S1 we derive transformations analogous to Pfaff's and Euler's for the func- 
tion 3^2(11,^2,6+ l;6i,e;x), involving nonlinear parametric constraints. Bailey 
and Slater (see, e.g., Slater [19, § 2.4.2]) found an Pfaff-like identity satisfied by 3F2, 
with three free parameters, but our identity, which extends theirs, has four. As a 
corollary, we extend other results of Bailey; e.g., his result that any very well 
poised 3F2 can be written in terms of a well poised 2-F1, and vice versa. 



2. Riemann P-symbols 

Any P-symbol tabulates the characteristic exponents of a linear (homogeneous) 
Fuchsian differential equation (FDE) Lu = on P . It partially characterizes the 
space of local solutions. The columns of the P-symbol list the exponents associated 
to each singular point of L. For instance, the GHE and HE have P-symbols 
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where the column and exponent orderings are not significant. Ordinary (non- 
singular) points may optionally be included, each with its own column. Any finite 
ordinary point of a second-order FDE has exponents 0, 1, and the point x = 00, 
if ordinary, has exponents 0,-1. Having 'ordinary' exponents is not sufficient for 
ordinariness, since zero-exponent Frobenius solutions may be logarithmic. 

Any second-order FDE on P 1 with n ^ 3 singular points is determined up to 
equivalence by its 2n exponents and n — 3 accessory parameters. Though the latter 
parameters are not shown, P-symbols facilitate the computation of transformed 
FDEs, as the following indicates. Automorphisms of P 1 are Mobius transforma- 
tions (also called homographies or linear fractional transformations) of the form 
M(x) = (Ax + B)/(Cx + D), where A,B,C,D e C and AD ^ BC. Under such 
transformations, singular points are accompanied by their exponents. That is, if 
M: P* -> Pi,, applying the change of variable x' = M(x) to an FDE Lu = on P*, 
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lifts it to (M*L)v = on P*, where if £ has singular points a', b', c' , . . . , with re- 
spective exponents ai, oli\ Pi, P2', 71, 72, ■ ■ ■ , the operator M*L will have P-symbol 
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Fuchs's relation states that the sum of the 2n exponents of any second-order FDE 
with n singular points equals n — 2. Mobius lifting trivially preserves this. 

Since there is a (unique) Mobius transformation that takes any three distinct 
points to any other three, one can normalize any FDE with n ^ 3 singular points by 
moving three of them to 0, 1, 00. Multiplying the dependent variable by (x — xo)~^ 
will decrement the exponents at x ~ xq by C, and those at x — 00 by — £; i.e., 
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Using ([5]), one can further normalize the FDE by shifting one exponent at each 
finite singular point to zero. For n = 3, 4, the resulting FDE is the GHE, resp. HE. 

Any nonconstant rational function x' = R(x) specifies a (deg i?)-sheeted cov- 
ering R : Pj. — > P^,. Applying the change of variables x 1 = R(x) to an FDE Lu = 
on Pi, lifts it to an FDE (R*L)v = on P*. 

Proposition 2.1. (i) The characteristic exponents of R*L at any point xq £ 
R~ 1 x' equal those of L at x' , multiplied by the multiplicity with which xq is mapped 
to x' . (ii) // the exponents of £ at x' are 0, l/k for some integer h > 1, and xq is 
mapped to x' with multiplicity k, then xq will be an ordinary point of R*L. 

Proof. Part (i) follows from the definition of a characteristic exponent. To 
prove (ii), it suffices to show that the zero-exponent Frobenius solution of (R*L)v — 
at x = xq cannot be logarithmic. By the definition of a lifting, this solution 
v = v(x) (normalization being arbitrary) can be written locally as u(R(x)), where 
u = u(x) is a zero-exponent Frobenius solution of Hu = 0. But such a solution can 
be logarithmic only if the other exponent is a positive integer; and here it is l/k. □ 

By examination, lifting L to R*L preserves Fuchs's relation. But under a 
non-M6bius rational lifting, singular points may disappear or appear. Part (ii) of 
the proposition indicates how an inverse image of a singular point may fail to be 
singular. Also, if x' is an ordinary point of £, with xq G R~ 1 x' mapped by R to x' 
with multiplicity k > 1, then x$ will be a singular point of with exponents 0, k. 

With the aid of Prop. 12.11 one can conjecture and verify rational transforma- 
tion laws for the special functions 2-F1 and HI. However, in any such identity the 
accessory parameter q of the HE will need to be considered separately. It should 
be mentioned that there are types of Hcun identity, such as differential transforma- 
tions, on which P-symbols throw less light. One example is the following. 

Theorem 2.2 ([20 ,§3). For every integer N ^ 0, one has the differential Heun 
identity 

D N Hl(a, q;l-N, P; 7, 8; x) oc Hl(a, q'; 1 + N, P + N; 7 + N,8 + N; x), 



where q' :=q + N(N - l)(a + 1) + N [(a + 1)7 + aS + e] 
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This result can be strengthened to an equality between HE P-symbols, namely 
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which says that any local solution of the left-hand HE, differentiated N times, 
is one of the right-hand HE. The P-symbol identity ^ is plausible, since the 
exponents on left and right formally correspond. For instance, at x = the left- 
hand Frobenius solutions are of the form x° and x 1-7 times analytic functions, each 
nonzero at x — 0; so (generically only!), their iV'th derivatives are of the form x° 
and x 1 ~(' 1+N ' > times such analytic functions. But a satisfactory proof of ^ along 
these lines would require a study of the series coefficients of the Frobenius solutions, 
since if any non-leading coefficient were zero, it would yield non-generic behavior. 
And of course, q' needs to be computed separately. In the following sections the 
focus is on non-differential identities, which are easier to conjecture and prove. 



3. Automorphisms of the GHE and HE 



P-symbols of normalized FDEs with n ^ 3 singular points can be written in a 
symmetrical way. Without loss of generality, the singular points may be taken to 
include x — 0,oo; and typically, x — 1 as well. Each finite singular point may be 
taken to have a zero exponent. The P-symbol of such an FDE will be of the form 
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where a\ — (and usually, 02 = 1). The canonical FDE solution F at x = a\ = 
will be written F(a,2, ■ ■ . , a n -V, Qj 0i, ■ ■ ■ , #n-i, #00; x), where q is a suitably defined 
vector of n— 3 accessory parameters. Necessarily a — (n—2 — 6\ — - ■ ■ — 6 n -i—8 QO )/2. 

The automorphism group of such normalized FDEs, or P-symbols, contains a 
subgroup of automorphisms that negate exponents, or more accurately exponent 
differences . These are called F-homotopies [18) . The subgroup is generated by 
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1, . . . , n — 1, and by the trivial a <-> (3 interchange at x — 00, i.e. 
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The group of F-homotopies is isomorphic to (Z2)™; or (Z2) Tl ~ 1 , if a «-> (3 is excluded. 

An F-homotopy at any of the n — 1 singular points other than x = a\ = yields 
a transformation of the canonical local solution F = F(x) at x = a\ = 0. Taking 
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the usual normalization F(x = <zi = 0) = 1 into account, one deduces 

(8) F(a 2 , ■ • ■ , a„„i; q\ 1 , . . . , 9, h . . . , n -i,0oc', %) 

= (1 - x/ai) 6i F(a 2 , ■ ■ ■ ,a n -x; q'; Ox, . . .,-6i, . . . , 6 n -\, #oo; x), 

i = 2, . . . ,n — 1, where for each i, the transformed vector q' of accessory parameters 
must be computed separately, by manipulating the FDE. One also has 

(9) F(a 2 , ■ • • , On-i; q; 6x, ■ ■ ■ , n -i> ^ooj x) 

= F(a 2 , ■ ■ ■ , a„_i; q; 6\,..., n -x, — B^; x). 

The group of F-homotopic transformations of F is isomorphic to (Wi 2 ) n ~ 1 ; or 
to (Z 2 ) n ~ 2 , if the trivial transformation ((9]), i.e., a <-»/?, is excluded. 

The GHE case n = 3, in which there are no accessory parameters, illustrates 
this. The identity {SJ, when n — 3 and i — 2, says that 

(10) 2 F!(«, /?; 7 ; x) = (1 - 2^(7 - a,j - 0; 7; 1), 

and © simply says that 2 -Fi(o! ! /3;7;a;) = %F\{0, a; 7; x). Equation (JTUJ) is Euler's 
transformation of 2 Fi [2], which is an involution that commutes with a <-> 0. 

The automorphism group of normalized FDEs Hu = with n singular points, 
or strictly of their P-symbols, also contains a subgroup of Mobius automorphisms, 
based on Mobius transformations. If Lu = is to be lifted to (M*L)v — along 
M: — > P^,, the inverse images of the singular points x' = 0,1, 00 of L must be 
chosen to be singular points of M*L. So there are n(n — l)(n — 2) possibilities, 
and there is no additional freedom in the choice of M. But since the order of the 
remaining n — 3 columns in the lifted P-symbol is not significant, one may permute 
them arbitrarily. So the Mobius subgroup is isomorphic to the symmetric group S n . 

The cases n = 3,4 illustrate this. In the n = 3 GHE case, x' = M(x) may be 
any of 3 • 2 • 1 = 6 homographies. The correspondence to permutations of 0, 1, 00 is 

x x 1 — x x ~ 1 1 1 

(o)(i)(oo), (of(:M, (oi)M, (oral), (o4(i), (01^). 

In the n = 4 HE case, x' — M (x) may be any of 4 • 3 ■ 2 = 24 homographies, given 
in [14] . The six that stabilize x = 0, which are of primary interest here, are 
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x 

(0)(la)(oo), 

1 

a ' 



x 
x — a 

(0)(laoo) 
1 

l-a ' 



Each M corresponds to a permutation of the HE singular points 0,1, a, 00, i.e., 
columns, as shown. But the permutation is to be interpreted in a special way. For 
instance, (0)(laoo) means 1 — s- 0, 1 1— > a', a 1— > 00, 00 1— > 1, for some a'. In general, 
aeP], and a' £ P*, differ. The value of a', in terms of a, is given in the third row. 

A Mobius automorphism based on a lifting along x' = M(x) must in general 
include an exponent shift, of the type shown in ([S]), in order to keep one exponent 
zero at each finite singular point. The formulas for the Mobius tranformations of 
the local solution F — F{x) at x — will include such factors. Any of the above 
n{n — l)(n — 2) Mobius transformations that fixes the point x = a\ = will yield 
a transformation of F. If the map x' = M(x) is not of the form x' — Ax, it must 
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be of the form M(x) = Ax/[x — M 1 (oo)], where M 1 (oo) is one of a 2 , . . . , a n -i- 
Taking F(x = ai = 0) = 1 into account, one obtains a transformation of the type 

(11) F(a 2 , ■ ■ - ,a„_i; q; 9\, . . . ,6 n -\,0oo\ x) 

= [1 - x/M-\oo)] ~° F(a' 2 , 4, . . . , a;_ i; </; ^, . . . , Ci, M(a;)), 

in which the transformed vector q' of accessory parameters must be computed 
separately. Any such formula is really (n — 3)! formally distinct transformations 
of F, since the n — 3 parameter pairs (03,03), ■ ■ ■ 1 ( a n-i,@n-i) on the left-hand 
side may be permuted arbitrarily. (Equivalently, the columns of the left-hand P- 
symbol, other than the x — 0,1, 00 columns, may be permuted.) The group of 
Mobius transformations of F is isomorphic to S n -±. 

The n = 3 GHE case is illustrative. Of the 3-2-1 = 6 homographies above, 
the only two that fix x = are the identity x 1— > x and x 1— > — they make up a 
group isomorphic to S 2 . By examination, the latter yields the identity 

(12) 2 F 1 (a,(3;r,x) = (l-x)- a 2 F 1 (a, 1 -(3; T , s£r), 

valid near a; = 0. This is Pfaff's transformation of 2 Fi [2], which is an involution. 

Any Mobius automorphism permutes the columns of the P-symbol, and any 
F-homotopic one acts on columns individually. So the S n subgroup of the former 
normalizes the (Z 2 ) n subgroup of the latter, and the full group of automorphisms 
is isomorphic to a semidirect product (Z 2 )" x S n . It is easy to see that this is the 
wreath product 1 2 \S n , the group of signed permutations of an n-set. Any element 
of Z 2 I S n = (Z2)' 1 xi S n is an ordered pair (h t a), comprising a Mobius and an 
F-homotopic automorphism, which act from right to left. The first permutes the 
singular points, i.e., columns of the P-symbol, according to a <E S n , and the second 
negates zero or more of the exponent differences 81, . . . , 9 n . Any element (h, a) can 
accordingly be written in a sign-annotated version of disjoint cycle notation. For 
instance, [0+] [l + a_oo_] signifies the signed permutation 1— > 0, 1 1— > a', a 1— > 00, 
00 1 ► 1 , followed by negations at x = 00 and x = 1. 'Even-signed' permutations, 
such as this one, include an even number of minus signs; they play a special role. 

Theorem 3.1. The automorphism group of normalized second-order FDEs on 
P 1 with n 3 singular points is the Coxeter group H n := Z 2 lS n , the group of signed 
permutations of an n-set. If the transformation a j3 is not explicitly included as 
an element, the automorphism group will be the Coxeter group T) n := \L 2 I S n ] , 
the group of even-signed permutations of an n-set. The transformation group of the 
canonical solution F — F(x) is Un-i, or 35 n -i if ot*-* f3 is not explicitly included. 

Remark. The orders of these four groups are |23 n | = 2 n n\, |D n | = 2" _1 rt!, 
|S„_i| = 2 n - 1 {n-l)l, and |D„_i| = 2 n " 2 (n-2)!. It is readily verified that £3,2)3 
are isomorphic to Z 2 x S4 and S4; and that H 2 , r D 2 are isomorphic to the dihedral 
group of order 8 and the Klein four-group (Z 2 ) 2 . The groups ® n ,D n are usually 
not defined as permutation groups, but as reflection groups consistent with certain 
Coxeter graphs [10] . Due to this, the notation used here is slightly nonstandard: 
our D 3 is usually denoted A3, and our T) 2 is not usually called a Coxeter group. 

Proof. This theorem follows from the preceding discussion, the transforma- 
tion group of F comprising all automorphisms that fix the singular point x = 0, 
and perform no F-homotopy there. If a <-> /?, i.e., Ooo 1— > —6*00, is not explic- 
itly included, then one may modify each automorphism with an odd number of 
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F-homotopies (negations) by following it by 9oc >— * —Qoo, to coerce even-signedness. 
This alteration is innocuous, since it leaves the FDE invariant. □ 

Example (n = 3). The correspondence between the elements of D 2l the group 
of even-signed permutations of the nonzero singular points x = 1, 00 of the GHE, 
and transformations of 2 F\ is 

(13a) [l+][oo+] ~ 2-F , i(«,/3; 7; x), 

(13b) [l_][oo_] ~ {i-xy-o-^F^-a^-frr,*), 
(13c) [l+oo+] ~ (l-x)- a 2 F 1 (a, 1 -0; 7; ^), 

(13d) [l_oo_] ~ (l-x)-^ 2 F 1 { 1 -a,p- 7; 

Here (| 13b>[) . (]1 3c|) are Euler's and Pfaff's transformations, and (jl3d[) is a twisted 
Pfaff 's transformation in which the first and second arguments, and also a and /3, 
are interchanged. Substituting any of (|13b|l . (fT3c|) . (|13dj) for F(a,b;c;x) is an invo- 
lutory operation. Any two of these three operations commute, and their product is 
the third, confirming that D 2 is isomorphic to the Klein four-group. If the trans- 
formation group of 2F1 is extended by explicitly including the involution a /3, 
the group becomes not D 2 , but !B 2 . The nontrivial coset of D 2 in % 2 comprises 

(14a) [l+][oo_] ~ 2 Fi{p,a\r,*), 

(14b) [l-][oo+] ~ {l-xy- a -P 2 F 1 ( 7 -/3, 7 7; x), 

(14c) [l+oo_] ~ (l-a;)- /3 2F 1 (/3,7~a; 7; ^), 

(14d) [l_cx>+] ~ (1 -x)- a 2 F 1 ( 1 -l3, a; 7; ^). 

Examination confirms that as expected, the set of transformations of 2 Fi(x) in 
(fl3)) .(|14 p . under composition, is isomorphic to the order-8 dihedral group. 

The automorphism group of the GHE is isomorphic to ©3 = S4 (of order 24), 
if a «-> (3 is not included. The group T) 2 yielding (I13ap - (|13d[) is a proper subgroup; 
a GHE automorphism may be more general than a transformation of 2 F\. In all, 
D3 yields 24 solutions of the GHE, which are the well-known solutions of Kummer. 
They are partitioned into 6 equivalence classes of size 4, one of which is p3a[) - (|13d[) . 

Example (n = 4). For any of the six homographies x' = M(x) given above, 
which in effect permute the nonzero singular points x — 1, a, 00 of the HE, one can 
rewrite (1111) in HE notation, obtaining a Mobius transformation of HI. The painful 
part is computing the transformed accessory parameter. The usual parametriza- 
tion (|lb[) of the HE is not adapted to this, and a closely related one is better. 

First, write q = a(3Q, where Q is an alternative accessory parameter. (Q, not q, 
was Heun's accessory parameter; the present convention was introduced by Erdelyi.) 
Second, perform a change of variable, x — a/x, i.e., (Ooo)(la), accompanied by an 
F-homotopy that restores the transformed HE to HE form. In P-symbol terms, 



(15) x° 



' 


1 


a 


00 


X 











a 




. 1-7 


1- s 


1 -e 


a 





' 


1 


a 


00 


a/x 











a 




(3 — a 


1 -e 


1 -S 


a — 7 + 1 
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That is, if F = v(x) is a solution of the HE (jlbj) . then x a v(x), regarded as a function 
of x — a/x, i.e., as a multiple of x~ a v(a/x), will satisfy a HE of the form 

a/3(x — Q) 



(16) D 



7 6 
x x — 1 



D s 



x(x — l){x - 



F = 0, 



with parameters a = a, /? = a — 7 + 1, j — a — (3 + 1, 5 — e, e = 5. By direct 
computation, Q = [0Q + (e - /3)a + (5 - /?)] /(a - 7 + 1). 

In computing the Mobius transformations of HI , 5 = a/x is a more natural 
independent variable than x, and Q a more natural parameter than Q or g. This 
is because in terms of x, each of the six homographies is affine linear, since it fixes 
the point X = 00, i.e., x — 0. From each homography x' = M(x), the afhne map 
x' = M(x) follows from x 1 = M(x) = a! jx' = a' /M(a/x). The correspondence is 

M (x) : x, t*, 2, -5- 

__ V / 5 X— 1 ' 2 — a ' a(:c — 1) ' a' :c — a 

A/ZYt^ ■ a-g 1 _ ^ a-S 2 1 — S 

iU V x ;- X J a_l ! IX, a ' a' 1-a' 

And from (|16|) . under any of these affine changes of variable i' = M(x), it is clearly 
the case that Q' = M(Q). So, Q' in terms of Q is quite easy to calculate. 

The homography M(x) = -zrp which comes from the permutation (loo) (a) of 
the nonzero singular points x = l,a,oo, is typical. For this homography, Eq. pip 
says that in a neighborhood of x = 0, Hl(a, q; a, 0;^, 5; x) can equivalently be 
written as (1 — x)~ a Hl(a' , q'; a', 0'; 7', 5'; ^zj), in which a' = a' = a, /3' = 
a — 5+1, 7' = 7, and <5' = a — /3 + 1, as is easily verified. What remains to be 
determined is the transformed parameter q' . But Q' — (a — Q') / (a — 1), since the 
affine linear map associated to M is M(x) — ^5f- Taking into account that 

Q = [0Q+{e- 0)a +(8- 0)) /(a ~ 7 + 1), 

Q' = [/3'Q' + (e' - 0')a' + (5' - 0')} /{a 1 - 7' + 1), 

and that q — a0Q and q' — a'0'Q', one finds q' = [—q + -faa)/(a — 1). 

The preceding procedure can be carried out for all six homographies. It yields 
the following bijection between (positively signed) permutations of the nonzero 
singular points x = 1, a, 00 of the HE, and Mobius transformations of HI. 

[l + ][a+][oo+] ~ Hl(a,q; a,0; j,S; x), 

[l+oo+][a+] ~ (l-x)- a Hl(^,^±^;a,a-5+l; 1 ,a-f] + l- ^), 

[l+][a+oo+] ~ (l-f)- Q ffi(l-a,-g + 7a;a,-/3 + 7 + ^7^; 

[l+o+][oo+] ~ FZ(i4;a,/3;7,a + /3-7-<5 + l; f), 

[l +a+0 o + ] ~ (i_|)-«^|( I l_ j ^;a > - J 9 + 7 + *;7,a-/3 + l; 5 § s ) ) 

[l + oo + a+] ~ (i-^-^^.^trsa; 

a,a-5+l;7,a + /3-7-5 + l; ^~^ ). 

Here, the second transformation is the just-derived one. It and the third are Pfaff- 
type transformations of Hl(a, g; a, 0; 7, (5; x). The fourth is a simple l 6 <-> e' one, but 
the fifth and sixth are quite novel, having no 2-F1 analogue. They cyclically permute 
the nonzero singular points, so they are not involutions: they are of order 3. 

The full transformation group of HI is isomorphic to ©3 = S4, or to S3 = 
Z 2 x S4 if the a <-> interchange is included. It is generated by the Mobius 
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transformations, and by F-homotopic ones of the type ([8]), which are easier to 
work out explicitly (details are left to the reader). The simplest such is 

[l_][a+][oo_] ~ (l-x) 1 ' 5 Hl(a,q-(S-l)ja; (3 - S +1, a - 5 +1; 7,2-5; x), 

coming from an F-homotopy at x = 1, which is analogous to (|13b[) . Euler's trans- 
formation of 2-Fi- An explicit bijection between the order-24 group D3 and the 
transformations of Hl(a,q;a, (3;^f,S;x) is given in [141 Table 2]. The 24 trans- 
formations were obtained by machine computation. But as the present treatment 
makes clear, they can also be obtained reliably by hand. 

The automorphism group of the HE is isomorphic to 2) 4 (of order 192), or to "84 
(of order 384), if a <-*■ (3 is included. The entire set of 192 HE solutions appears in 
Ref. [14] . and the bijection to D4 is given explictly. The 192 are partitioned into 
8 equivalence classes, one of which comprises the just-mentioned 24 transformed 
versions of HI. The remaining 7 come from the cosets of D3 in D4. For later use, 
one of the 168 = 192 — 24 solutions of the latter sort will be given here; namely, 

(17) [a + 0+l + ][cx) + ] ~ Hl(^,^^;a,(3;a + P-y-5 + l, r ,^). 

This is obtained from a Mobius lifting M : P* -> P*,, where M{x) := ^p. Since 
M" 1 (0) = a 7^ 0, it is not a transformed version of HI: it is a zero-exponent solution 
of the HE near the singular point x = a. It and Hl(a, q; a, /3; 7, 5; x) do not equal 
each other on their common domain of definition, if this domain is nonempty 



4. Quadratic and quartic transformations of HI 



The transformations of HI derived in §|3]were based on Mobius transformations 
of P 1 , i.e., on degree-1 covering maps M: P* — > But if R: — > P*, is any 
rational map, a HE Lu = on the Riemann sphere P*, will lift along R to an 
FDE (R*£j)u — on P*. The characteristic exponents of the operator R*L are 
determined by Prop. [2~T1 If R*L has only four singular points, the lifted FDE will 
be a HE in its own right. Such a lifting is said to be a degree- (deg R) transformation 
of the HE. If, moreover, R(0) = 0, then a formula relating a local solution (at x = 0) 
of the form Hl(x), to one of the form Hl(x') — Hl(R(x)), will result. 

Erdelyi [5J §15.3] pointed out the existence of quadratic and quartic transforma- 
tions of the HE, though he did not derive the corresponding functional equations 
satisfied by the special function HI. He showed that if any two of the singular 
points x' = 0, 1, a', 00 of the HE on PjJ., have exponent difference 5, then there is a 
quadratic lifting to a HE on P* . (The singular points x = 0, 1, a, 00 of the lifted HE 
may differ, since a' ^ a in general.) For instance, one has 



(18) (1-xY 



' 


1 


a 


00 


x 











2a 




. 1-7 


7 — 2a 


1-7 


7 





in which the map x' — R(x) is of the form 






1 


a' 00 


R(x) | 








a 




1-7 


1 

2 


1 7-« 





(19) 



R(x) = A 



x (a — x) 
1-x 
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provided that a, a' 6 P 1 \ {0, 1, oo} are related by 

(20) a 2 (l-a') 2 -16(l-a)a' = 0, 
and the 'multiplier' A — A(a, a') is given by the formula 

(21) A = A(a,a') L + " 



2(2 - a) 



The P-symbol identity (HHJ) follows from Prop. [HEJ since J? _1 (0) = {0, a} and 
i? _1 (oo) = {l,oo}; and as one can check, x' = l,a' are critical values of the map 
x' = R(x), the corresponding critical points being mapped doubly to them. The 
latter requirement, which causes the singular points x' — l,a' to 'disappear' when 
lifted to pi, is the source of the quadratic constraint (|2"0|) and the formula (|21|) . 

The Riemann surface defined by (|20[) is of genus zero. Equivalently, a, a' can 
be expressed in terms of a complex parameter t. By examination, one can choose 

(22) ait )= t l±R, / (t) = _^ ; A (t)= 



(t + 4) 2 ' w (i + 8) 2 ' w v * + 

The choices t = 0, —4, —8, oo are 'unphysical', since each corresponds to one of a, a' 
equaling one of 0, l,oo. So, the possible quadratic Heun liftings of the form (fT5|> 
are parametrized by the quadruply punctured sphere Pj \ {0, —4, —8, oo}. 

The quadratic (and biquadratic quartic) transformations of HI can now be 
determined. To avoid degenerate cases, it will be assumed that the underlying 
rational map R: P* — > P*, has the property that each critical value, on P*,, is one 
of the HE singular points x' — 0, 1, a, oo. 

Theorem 4.1. There is a unique quadratic transformation of HI, up to pre- 
and post-composition with Mobius and F-homotopic transformations of HI. It is 

Hl(a, q; 2a, 7; 7, 2a — 7 + 1; x) 

= (1 - x)~ a Hl(a' , A{q - jaa); a, 7 - a; 7, \ ; R(x)), 

where a, a' are related by (|20[) , and the quadratic map R and multiplier A are defined 
by (|19[) . (|21[) . This equality holds on a neighborhood of x = 0. 

Proof. Any quadratic map R must have branching schema 1 + 1 = 1 + 1 = 2 = 
2. (A schema of this sort, with four 'slots', lists the multiplicities with which the 
inverse images of the four singular points x = 0, 1, a', 00 on P*, are mapped to them; 
orderings of slots and inverse images are not significant.) In general, the lifted HE 
will have six singular points. The number will be reduced to four, by Prop. 12.11 if 
each of the two singular points on Pj., listed last has exponent difference |. 

Let the schema be annotated by underscoring each inverse image that is a 
singular point of the operator R*L on P*. Then there is a unique annotated schema 
corresponding to a lifting of an HE to a HE, namely 1+1 = 1+1 =2 = 2. For this to 
yield a transformation of HI, one must have R(0) =0. So x' — must be assigned 
to one of the first two slots. The remaining singular points x' = l,a',oo can be 
assigned arbitrarily to the remaining three slots; and the singular points x = 1, a, 00 
to the remaining inverse images of type 1. Any such arrangment can be reduced by 
pre- and post-composition with Mobius transformations to the scheme (|18|) , except 
that F-homotopic transformations may also be applied. 
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The HI identity of the theorem can be read off from (fTSj) , except for the value 
q' = A(q — ^aa) of the right-hand accessory parameter, which follows by direct 
computation: applying the change of variable x 1 = R{x) to the HE Lu = 0. □ 

Remark 4.1.1. The full set of quadratic transformations of HI is obtained 
by (i) composing with Mobius maps, i.e., using M2 o Ro Mi rather than R as the 
covering map (the left side of the identity being transformed by M\~ , and the right 
by M2), and (ii) acting on either side by F-homotopies. Quadratic transformations 
are therefore bijective to the elements of D3 x CD3, modulo an appropriate stabilizing 
subgroup. (Here 'a «-> /?' is disregarded.) A complete list will appear elsewhere. 

Remark 4.1.2. To place Theorem l4.1l in context, there is no analogous unique- 
ness result for the quadratic transformations of the Gauss hypergeometric func- 
tion 2-Fi- There is a unique annotated schema for a quadratic lifting of a GHE 
on P*, to a GHE on P*; namely, 1 + 1 = 2 = 2. But there are two ways of as- 
signing x' — 0: it can be assigned to the first slot or the second, and the resulting 
quadratic maps R(x) are distinct (up to pre- and post-composition with Mobius 
transformations, i.e., with Pfaff's involution x t— ► x/(x — 1)). The two choices yield 

ArrfA — <r\ ~ 4x 4x 4x(l-x) . 

nyx) — ^l 1 ^ x), (1-xJ 2 ' (1+x) 2 ' (l-2x) 2 ' 

R ( X ) = (2-x) 2 ' 4(l-x) ' 

respectively. Taking pre- and post-composition with F-homotopies (i.e., with Eu- 
ler's transformation) into account, one sees that the transformations coming from 
the first alternative are bijective to the elements of the group D2 x ^2, and those 
from the second are bijective to the same group, modulo a stabilizing subgroup of 
order 2. A partial tabic of quadratic transformations of 2-F1 is given in [HI §2.11]. 
The existence of two distinct types has been stressed by Askey [3]. 

Erdelyi also pointed out that if three of the exponent-differences of the HE 
on V\i are equal to h, such as those at x' = 1, a', 00, then it can be lifted to a HE 
along three distinct degree-2 covering maps x' = R(x), and each can be followed 
by a second degree-2 lifting, leading to a biquadratic quartic transformation. By 
examination, any such composition of two quadratic maps, x' — S(x), has branching 
schema l_+l_+l_+l_=2 + 2 = 2 + 2 = 2 + 2. This schema is exemplified by 

4ax(l - x)(a - x) (a-2ax + x 2 ) 2 a(a - 2x + x 2 ) 2 

( 3) S[X}= (a~x 2 ) 2 =1 (a-x 2 ) 2 =a (a~x 2 ) 2 ' 

for which a = a'; for this map, 5 _1 (0) = {0, 1, a, 00}. The schema is so symmetrical 
that when deriving quartic transformations of HI from it, there is essentially only 
one way of assigning x' = 0, 1, a', 00 to the slots, and x — 0, 1, a, 00 to the inverse 
images of type 1. Proceeding as above, and computing the accessory parameter q 
of the lifted HE from q 1 by direct computation, one obtains the following. 

Theorem 4.2. There is a unique biquadratic quartic transformation of HI, 
up to pre- and post-composition with Mobius and F-homotopic transformations 
of HI. Valid on a neighborhood of x — 0, it is 

Hl(a,q: 27- 1,7; 7,7; x) 

= (1 - x 2 /a)-^ Hl(a, q/4; 1 - |, 2 + I; 7, |; S(x)), 

where the biquadratic quartic map S is defined by (|23D . 
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If 7 = 1/2, then all four of the exponent-differences of the HE are equal to h. 
In this quite degenerate case, the biquadratic HI identity specializes as follows. 

Corollary 4.2.1. The special function H(a.q;x) := Hl(a, q; 0, h; \, ^; x), 
where a S C \ {0, 1} and 9 £ C, which is defined on \x\ < min(l, |a|), satisfies 
the functional equation H(a, g; x) = H(a, g/4; S(x)) on a neighborhood of x = 0. 

Remark. By clever changes of variable, due to L. Carlitz (see Valent [211 ). one 

where sn _1 is the inverse of the 



can show H(a,q;x) = cosh 2^J q/a sn 1 (y / ir) 
Jacobi elliptic function sn(- 1 1/a), and the branch that vanishes at x = is meant. 
So the corollary says that sn _1 (y / x) = sn -1 S(x)) /2, i.e., that sn 2 (2u| 1/a) = 
5" (sn 2 (u 1 1/a)). That is, the 7 = 1/2 case of the biquadratic HI transformation is 
really the duplication formula for the Jacobi function [§] §13.17], in disguise. That 
the biquadratic quartic function S plays a dual role is not well known. 

5. From HI to 3 F 2 

The local Heun function Hl{x) := ^2^L c(n)x n has been less extensively in- 
vestigated than the Gauss hypergeometric function 2Fi{x), or even the generalized 
hypergeometric functions %F2(x), ^F^lx), etc. The presence of a fourth singular 
point x — a in the HE (|lb|). or equivalently, the fact that the recurrence (|2bj) sat- 
isfied by the coefficient sequence c: N — > C is of second order rather than first, is a 
novel and somewhat disconcerting feature. It is natural to ask whether, for certain 
parameter choices, HI can be reduced to more familiar transcendental functions. 

That HI may reduce to 3F2 was shown by Letessier et al. In this section 
their result is rederived from scratch, starting from the l Hl side' rather than the 
'3F2 side,' as it were. A connection to P-symbols will be made. As usual, it will be 
assumed that the parameter 7 is not a nonpositive integer. 

First, consider how the singular point x = a of (jlb[) can be removed or otherwise 
tamed. The simplest case is e = 0, when the characteristic exponents of x = a are 
0,l — e=l, which are those of an ordinary point. (The case e = 2, which is related 
to this one by an F-homotopy, is left to the reader.) Since the difference between 
the two exponents is an integer, the local solution corresponding to the smaller one, 
zero, will generically be logarithmic. But a glance at (jlb[) reveals that if q — a(3a 
when e = 0, then the HE will reduce to the GHE. Since 5 = a + j3 — 7 — e+1, one 
has the reduction Hl(a, a/3a; a, /3;7, a + f3 — 7 + 1; x) — 2-Fi(a, f3; 7; x). 

The next simplest case is e = — 1, when the exponents of x — a are 0, 2. (The 
case e = 3, related to this one by an F-homotopy, is also left to the reader.) As 
before, the zero-exponent local solution of the HE at x = a will generically be 
logarithmic. If it is not, then by definition x = a will be a (simple) apparent 
singular point: all local solutions of the HE at x = a will be analytic. 

Theorem 5.1. At fixed a, (3, 7 S C, the family of HEs with e = — 1 (and hence 
5 — a + P — 7 + 2), the singular point at x = a being required to be apparent, can 
be parametrized by an auxiliary parameter e G C according to 

e(e- 7 + l) (e + l)(e-7 + l) 

■ — ap ■ 



(e ~ a)(e ~ (3) ' (e — a){e — j3) 

Proof. The singular point x — a will be apparent iff q satisfies a certain 
algebraic condition. It can be worked out by substituting F(x) — X)^Lo c(n)(x—a) n 
into (fib)) , with c(0) = 1, and solving for c(l), c(2), etc. In calculating c(2) a division 
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by zero will occur, unless q satisfies the condition, in which case the equation for c(2) 
will be • c(2) = 0. The condition is that the right side of this equation equal zero. 

More efficiently, one may exploit the explicit formula (| 1 T[) for the zero-exponent 
local solution of the HE at x = a. When e = — f , it reduces to 

(24) Hl(a', q'; a', 0'; 7 ', 6'; x') := Hl{*±, a, /?; -1, 7 ; 2=£). 

Since 7' = — f, this HI has exponents 0,2 at x' = = 0, and is generically 
logarithmic there. The condition on q' for Hl(a',q';a',(3';—l,5';x') to be non- 
logarithmic follows readily from the recurrence (|2b[) . One hnds that the equation 
determining c(2) is of the form • c(2) = 0, and therefore has a solution, iff 

(25) q' 2 + [(a' + /3'-5' + l) + (6' - l)a'} q' + a' [3' a' = 0. 
Substituting the expressions for a', q'; a', (3'; 5' from (j2"4")l into (|23|) yields 

(26) q 2 + [(7 - 1) - (2a/? + a + 0)a] q + a/3a [(a/3 + a + (3 + l)a - 7] = 0, 

which is the desired condition on q. 

If a, (3, 7 are fixed, Eq. (|26|) defines a curve in the a-g plane. (In the complex 
domain, it is a two-sheeted Riemann surface over .) Being quadratic, it is of genus 
zero and can be uniformized by a P^valued parameter, say e. By examination, the 
parametrizations a = a(e), q — q(e) supplied in the theorem will work. □ 

Remark 5.1.1. The theorem says that at fixed a,/3,7, the family of HEs for 
which x = a has exponents 0, 2, but is a relatively 'tame' apparent singular point, 
is parametrized by the Riemann sphere P*. Generically this sphere is sextuply 
punctured, since e = 0, 7 — 1 yield a = 0, e = a/3 /(a + (3 — 7 + 1), 00 yield a = 1, 
and e = a, (3 yield a = 00; and necessarily a 7^ 0, 1, 00. 

Definition 5.2. The special function G(x) :— G(a, (3; 7; e; x) is defined to 
equal Hl(a, q; a, (3; 7, 5; x) on a neighborhood of x = 0. Here (5 = a + /? — 7 + 2 
(so that e = —1), and a, g are given in terms of e (and a, /3, 7) in Theorem 15. II 

How the parametrized local functions G(x), 3^2(21) arc related will now be ex- 
plained. The following facts about 3-F2 (x) will be needed. It can be defined as a 
solution of a canonical third-order differential equation, or as a sum X)^o d(n)x n , 
where d: N — * C satisfies a first-order recurrence [H Ch. 4]. The differential equa- 
tion is [xD rii=i [pD + bi — 1) — x Y\^ =1 (xD + a,)] F — 0; or equivalently, 

(27) ( j3 + U±^j±^±3 _ bi + 62+l l D 2 
I x — 1 2(2; — 1) 

aia 2 + a 2 «3 + a 3 ai + ai + a 2 + a 3 + 1 bxb 2 

x(x — 1) a; 2 (x- — 1) 

Here a\, 0,3) b\, 62 are complex- valued parameters, and the equation is invariant 
under separate permutations of ai,a2,a 3 and 61,62- The singular points are x = 
0,l,oo, with respective characteristic exponents 0, 1 — b%, 1 — 62; 0, 1, s; 01, a,2, (13, 
where s := 61 + 62 — o-i — «2 — 03 is the 'parametric excess.' If either of 61, 62 is a 
nonpositive integer, the local solution of ([27]) corresponding to the zero exponent 
at x = will generically be logarithmic; but if neither is a nonpositive integer, 
which will be assumed henceforth, the solution will be analytic on \x\ < 1. When 
normalized to unity at x = 0, it is denoted 3^2(2;) := 3^2(01, 02, 0.3; 61, 62; x). 



D 



aia 2 a 3 
x 2 (x — 1) 



F = 0. 
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Equation (|27|) is not the most general third-order (linear, homogeneous) FDE 
with the above exponents, but it is the most general one associated to a first-order 
recurrence. The coefficient sequence d: N — > C of 3^2(3?) satisfies 



(28) 



{(n + bi)(n + b 2 )(n + 1)E - (n + a x ){n + o 2 )(n + 03)} d = 0, 



initialized by d(0) = 1. Here i? is the shift operator. This hypergeometric recur- 
rence is of course the traditional way of defining 3^2(2;)- 

Theorem 5.3. The local function G(a, (3\ 7; e; x) can be expressed in terms 
0/3-P2, as 3 F 2 (a,P,e+ l;j,e;x). 

Proof. It is readily verified that if e = — 1 (so that 5 = a + (3 — 7 + 2), and 
a, q are given the values supplied in Theorem 15. 11 then the second-order difference 
operator in the recurrence equation (|2b[) . the normalized solution of which is the 
coefficient sequence c: N — > C in the series expansion of G, has the factorization 



(29) -±— { < e ~r<+\ E - l\ 

V / n+e+1 1 (e— a)(e— p) J 



|(n + 7)(n + e)(n + l)^ - (n + a)(n + j3)(n + e + 



The right factor is the operator appearing in (|28|) . with cii, 02, 03; 61, 62 set to 
a, /?, e + 1; 7, e. So the hypergeometric sequence 0?: N — > C with these parameters 
is a solution of (|2bj) . Its initial conditions are those of c, i.e., d(0) = l,d(— 1) = 0. 
So c = d, and G{x) = 3 F 2 {x). □ 

Alternative Proof. It is readily verified that if at, <z 2 , 03; 61, b% are set to 
a,P,e + l;7,e, then the third-order differential operator in Eq. (|27|) . the unique 
(normalized) analytic local solution of which is 3-F2, has the factorization 



(30) 



D 



1 



1 



1 



D 



x-1 



D 



a/3 x — q 



x(x — l)(x — a) 



in which e:= —1, 5 := a + P — 7 + 2, and a, q are defined as in Theorem 15.11 The 
right factor is the operator appearing in the HE (jlb[) . So, 3^2 (%) — G(x). □ 

Remark. The operator factorization (|2"9"|) , resp. (|3U|) . can be viewed as taking 
place in the noncommutative polynomial ring generated by n, E, resp. by x, D. (The 
noncommutativity is specified by [E,n] — E, resp. [D,x] — 1.) Such factorizations, 
when they exist, can be found by noncommutative Grobner basis techniques [15] . 
There are also specialized algorithms for finding hypergeometric-sequence solu- 
tions to recurrences, i.e., for finding first-order right factors of difference opera- 
tors [6l 1 16] . But in general, it is difficult to characterize when a difference or dif- 
ferential operator will factor. If e = — 1 and (a,q) lies on the algebraic curve (|26|) . 
then the second-order difference operator associated with HI will factor, as shown 
in ([29]) . Conversely, the factoring of the third-order differential operator associated 
with 3^2(01 j &2, e + 1; b\, e; x), as in (|30p . yields a parametrization of the a-q curve. 

Corollary 5.3.1. G(a, e; x) can be expressed in terms of the Gauss func- 
tion 2-F1, as 2F\{a, [3; 7; x) + (a/3/je)x 2-Fi(a + l,f3 + 1; 7 + 1; x), or equivalently, 
in terms of 2-F1 and its derivative as e~ x [xD + e] 2^1(0;, /3; 7; x). It can also be 
written as [(e — 7 + l)/e] 2^i(a, 0; 7; x) + [(7 — l)/e] 2-Fi(a, /?; 7 — 1; x). 
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Proof. For any generalized hypergeometric function, the case when an upper 
parameter exceeds a lower one by an integer is special: it can be reduced. For 
instance (131 §5.2], 3-^2(01, fi2,e + l;bi,e;x) can be written in terms of the Gauss 
function, as 2-^1(01, ^2; b\\ x) + (aia2/&ie)a; 2-^1(01 + 1,«2 + + Ij^), i- e -, as 
e~ 1 [xD + e] 2^1(01, 0-2] br, x). This can be verified by comparing series expansions. 
It is also the case that the general function F := 3^2(^1, 02, 0,3; &i, &2! x) satisfies 
the three-term contiguity relation 

(31) F(a 3 +) = [(a 3 - b x + l)/a 3 ] F + [(h - l)/o 3 ] 

where 03+ and 61— indicate incrementing and decrementing of parameters [IT] . 
Setting (a 1, 0,2, 0,3; b±, 62) = (a, /?, e; 7, e) yields the second representation. □ 

By Theorem 15.31 the function G(x) has both an HI and a 3F2 representation. 
Their relationship is revealed by their P-symbols, which are respectively 
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00 
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00 
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- 7 


1-6 
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— e 


1 


1 


e+1 





with (5:=a + /3 — 7 + 2. In the third-order P-symbol on the right, the ordinary point 
x = a is listed explicitly, with exponents 0, 1, 2. (Any order-fc FDE has exponents 
{0, . . . , k — 1} at any finite ordinary point.) The final exponent row comes from the 
left-multiplication by a first-order differential operator that was performed in (|30|) . 
but the two P-symbols are otherwise the same. Appending the '1' to the x = a 
column converts x = a to an ordinary point, leaving only three singular points. 

This reduction of HI to 3-F2 is therefore an example of desingularization: re- 
moving an apparent singular point from a differential operator by left-multiplying 
by an auxiliary, appropriately chosen operator. Such desingularizations may be 
constructed algorithmically, as Abramov et al. [T| show. They also show how to 
desingularize difference operators, which in a certain sense have singular points 
too. Any order-d operator L := '^2'j = aPj{n)E : ' , where the pj are polynomials with 
Pd =/= 0, may be viewed as acting on functions defined on C n , the entire complex 
n-plane. The difference equation LF = will have d independent solutions that 
are meromorphic on C ra and analytic on a common left half-plane. The (leading) 
singular points of L on C„ are the roots of the polynomial pd{n — d), by definition; 
they are possible pole locations. If n = n is a root, but none of the d solutions 
has a pole there, it is said to be apparent. A desingularized version of L can be 
obtained by, e.g., left-multiplying L by an appropriate difference operator. 

The right factor in (|29p . the first-order difference operator associated to the 
function 3-F2(a, 0, e + 1; 7, e; x), has singular points at n = 0, 1 — 7, 1 — e. It can be 
shown that the one at n = 1 — e is merely apparent. The left-multiplication in (|29|) 
removes this apparent singular point. So, the reduction of this 3-F2 to HI is also a 
desingularization. The two desingularizations are dual to each other. 



6. 3F2 transformations 

The generalized hypergeometric function 3F2 :— 3^2(011, 02, 03; b\, 62; x) satis- 
fies many identities, but its theory is less developed than that of 2-Fi- It was shown 
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in §[5] that the ad hoc special function 

G(a, (3; 7; e; x) = 3 F 2 (a, 0,e+l; 7, e; x) 

— e -1 [x D + e] 2 -Pi (a, /?; 7; cc) 

has a Heun representation. As a result, certain transformation laws for G (and 
hence for a restricted 3F2) can be deduced from those for HI. 

Pfaff and Eulcr-likc transformations of this 3-F2 are given in Theorem 16.11 be- 
low, which is the main result of this final section. Each transformation has four 
free parameters and a nonlinear parametric constraint. Part ([1} of the theorem 
subsumes several interesting identities discovered by Bailey and Slater, which have 
fewer free parameters. They appear as corollaries of the theorem. 

Some terminology from general hypergeometric summation theory will be used. 
The function r+ iF r (x) :— r+ \F r (ai, . . . , a r+ i; b±, . . . , b r ; x) can be defined as the 
solution of a canonical order-(r + 1) FDE, or as a sum J2^Lo d(n)x n that converges 
on \x\ < 1. The FDE is [xD]J r %=1 {xD + b, - 1) - xU^(xD + ai )]F = 0, and 
d: N -> C satisfies [(n + 1) \X=i( n + h)E - Ui=l(n + a % )] d = 0. 

In general, if {a^ij 1 and {6;}[ =1 can be separately permuted so that a\ + 1 = 
d2 + i>i = ■ ■ ■ = a r+ i+b r , then the function r +\F r is said to be well poised. If each of 
these r + 1 equalities holds except the last, it is nearly well poised, or simply nearly 
poised. If r +\F r is well poised, resp. nearly poised, and also one of the {bi}^ =1 , 
resp. {piYiZx ) equals ai/2, it is very well poised, resp. nearly very well poised. 

Theorem 6.1. 

(1) The function 3F2 satisfies the Pfaff-like identity 

3F 2 (a 1 ,a 2 ,e+ 1; 61, e; x) 

= {l-x)- a > 3 F 2 ( ai , br -oa-l.e'+l; 61, e'; ^), 
m which the e <-> e' correspondence is 

, (61 - a 2 - l)e 
e = . 

e - a 2 

(2) TTie function %F 2 satisfies the Euler-like identity 

3 F 2 (ai, a 2 , e + 1; b\, e; x) 

= (1 - a;)*!-^-^-! 3^3(6! _ ai - 1, 61 - as - 1, e" + 1; &i, e"; i), 
in which the e <-> e" correspondence is 

„ (61 - ai - l)(6i - a 2 - l)e 

6 = . 

(&i — ai - a 2 - l)e + aia 2 

Here e, e', resp. e, e", are taken to be finite. Also, the Mobius transformation e 1— > e', 
resp. e h e", which is of the lower triangular form e 1— > gferg , is assumed to be 
nonsingular, which will be the case if all upper parameters are nonzero, i.e., if 
neither $F 2 is identically equal to unity. 

Proof. The transformations of Hl(x) := Hl(a,q;a,(3;^,8;x) analogous to 
Pfaff 's and Euler's transformations of 2 -F\ are, from §[3J 

[l+oo+][o+] ~ (l-x)- a Hl(^ 1 ,^±^-,a,a~-S+l- 1 ,a~f3+l-^- 1 ), 

[l_][a+][oo_] ~ (1 -xf- s Hl(a,q- {5 - l) 7 a; - 5+ I, a -5+ 1; 7, 2 - 5; x). 
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Theorem 15.11 gives the parametrization (a, q) = (a(e),g(e)) of the case when the 
singular point x = a is a (simple) apparent one. In this case, Hl(x) equals 
G(a, /3; 7; e; x) = ^2(0., /3, e + 1; 7, e; x). The transformed versions of Hl{x) have 
similar parametrizations (a',q') = (a'(e'), q'(e')) and (a",q") = (a"(e"),q"(e")) , 
and may likewise be written in terms of 3F2. Here a' — -^j, q' = ~ g +T Q!a ; and 
a" = a, q" = q — (5 — l)"/a. The transformations therefore yield isomorphisms be- 
tween genus-zero algebraic curves, (a,q) <-> (a',q') and (a,q) «-> (a",q"), which on 
the parametric level are performed by Mobius correspondences e «-» e' and e <-> e". 
Computing these explicitly yields the formulas in the theorem. □ 

Remark. The transformations of the theorem can also be verified by expressing 
each 3F2 in terms of two 2-Pi's, as in Corollary 1 5. 3. 1[ and exploiting results on 2-F1. 

The Pfaff and Euler-like transformations generate an order-4 group of transfor- 
mations of the family of 3^2 's of the form 3-F 2 (ai, a 2 , e + 1; 61, e; x). The group is 
isomorphic to ©2 = Z2 x Z2. (The third nontrivial element of this group is a twisted 
Pfaff transformation, obtained from the Pfaff-like one by interchanging <2i,a 2 .) If 
0-1 <-> a2 is included explicitly, the resulting group will be isomorphic to ^2, i-e., to 
the dihedral group of order 8. This group structure, familiar from 2-F1, is inherited 
from the automorphism group of the HE. 

COROLLARY 6.1.1. The function 3F2 satisfies the Pfaff-like identity 

3 F 2 (ai, a 2 , -a 2 + 1; h, -a 2 ; x) 

= (l- x )-^ 3 F 2 (a 1 ,b 1 -a2-l, b -^f^; h,^f^; ^). 

PROOF. Set e = -a 2 in Theorem ETMrj. □ 

Remark. This special case, with three free parameters and no nonlinear con- 
straint, was discovered by Bailey and Slater [19| (2.4.2.8)] and rediscovered by 
Baweja [4]. Their derivations were altogether different. The right-hand 3F2 is 
nearly very well poised, as one sees by permuting parameters. 

Corollary 6.1.2. There is a three-parameter reduction of 3F2 to a single 2-F1; 
namely, 

3 F 2 (ai,a 2 ,e + 1; 61, e; x) = (1 - x) 2 i<i(ai + l,a 2 + 1; 61; x), 

in which e :— e(oi, &2; 61) = 0,10.2/(0.1 + 02 — bi + 1). 

Proof. If e = aia^/(ai + 02—61 + 1) then e' = ai, reducing the right-hand 
side in Theorem WTH^i to (1 - x)- ai 2 -Fi(ai + 1, h - a 2 - 1; h; To obtain the 

given right-hand side, use Eq. (fl2|) . Pfaff 's transformation of 2-F1 ■ □ 

COROLLARY 6.1.2.1. Any very well poised 3 F2 can be reduced to a well poised 2 Fi, 
according to 

3 F 2 (a,f3,a/2 + 1; a- f3 + l,a/2; x) = (1 - x) 2 Fi(a + 1; a-0+l;x). 

PROOF. Set (01.09; 61) = (a,B:a - 3+ I) in Corollary 16". 1.21 □ 

Remark. This reduction of a very well poised 3F2 was also found by Bailey. 
It is sometimes reproduced in a different but equivalent form |19l (2.4.2.10)], in 
which Eq. (fT0|) , Euler's transformation of 2-F1, has been applied to the right side. 
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(l-s)a- (l+s)/3- (|-t) 



a- (3 



Corollary 6.1.2.2. For any s,t £ C, £/ie parametrized function 
(l-x) 20 - 1 3 F 2 (2a- 1, (l-s)a-(l+s)/3-(i-t), e+1; (l-s)a+(l-s)/3-(i-*), e; x) , 
zn which 

e := (a - ±) 

zs stable under the interchange a «-> /?. 

Proof. By Corollary 16. 1.2[ this function reduces to 

(1 - a;) 2Q a-Fi (2a, (1 - s)a - (1 + s)/3 + (± + *); (1 - s)a + (1 - s)j9 + (\ + 1); xY 

Applying Eq. (flT))) , Euler's transformation of 2-F1, has the effect of interchanging a 
and /3 in this expression. □ 

Remark. For each s,t £ C, the invariance under a <-> p is really a special case 
of Thcorcm l6.1t[2|) . the general Euler-like transformation of 3^2(01, 02, e+1; 61, e; x). 

Corollary 6.1.2.3. There is an involutory transformation of the two-parameter 
family of very well poised 3F2 's, namely the a <-> (3 interchange 



(l-x) 2Q - 1 3 F 2 (2a 



1, a 



2 ' 2 
(l-x) 2 ' 3 - 1 3^ 2 (2/3-l,/3 



a + ^ 
- a - 



2 3 2 3 xj 
i 13+ |; + a + i, /? - i; x). 



PROOF. Set s = t = in Corollary [6"X2~2l 



□ 



Remark. This transformation law was also discovered by Bailey (see [8j §4.5]). 
It can be derived directly from Corollary 16. 1.2.11 which reduces both sides to the 
same well poised 2-F1 • The generalization given in Corollarv l6.1.2.2l is new, as is its 
interpretation as a transformation of Euler type, a special case of Theorem I6.1lj 2j) . 

It should be noted that the two sides of this transformation formula are so- 
lutions of different third-order differential equations. The solution spaces of these 
equations have nonempty intersection, but are distinct. This phenomenon is clari- 
fied by the respective P-symbols, which by examination are 
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The first two rows of exponents are the same in both, but the third row differs. 
This partial agreement is what one expects to see in two third-order operators that 
are distinct left multiples of a common second-order operator. 
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